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in  analogy  with  the  O'Neil,  Winfrey,  Malmberg  (OWM)  problem  but£^wes  a spatially 
averaged  description  based  on  the  exact  energy  and  momentum  conservation  laws 
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relation.  This  model  shows  that  the  beam  can  be  clamped  in  velocity  with  the 
momentum  push  being  transferred  to  the  waves.  The  model  has  been  extended  to  the 
relativistic  and  multi-mode  cases.  For  E = 0 the  nonrelativistic  results  are 
in  good  agreement  with  the  work  oP\(OWM) , while  in  the  relativistic  case  the 
model  reproduces  the  nontrivial  features  of  the  computer  study  by  Lampe  and 
Sprangle.  \ 


ABSTRACT 


An  idealized  problem  is  investigated  which  illustrates  the  role  of  wave- 
particle  interactions  in  the  evolution  of  runaway  beams.  The  model  considers 
the  interaction  between  a weak  cold  beam,  driven  by  an  external  static 
electric  field  E,  and  waves  quantized  by  the  geometry.  The  waves  may  corres- 
pond to  Gould-Trivelpiece  modes  fixed  by  the  length  of  the  experimental  device, 
or  to  the  finite  Fourier  modes  encountered  in  computer  simulations.  The 
physics  consists  of  the  sweeping  of  the  accelerated  beam  through  the  resonance 
provided  by  each  cavity  mode.  This  process  is  formulated  in  analogy  with  the 
O'Neil,  Winfrey,  Malmberg  (OWM)  problem  but  uses  a spatially  averaged  descrip- 
tion based  on  the  exact  energy  and  momentum  conservation  laws  with  the  dynamics 
simplified  through  a WKB  representation  of  the  dispersion  relation.  This  model 
shows  that  the  beam  can  be  clamped  in  velocity  with  the  momentum  push  being 
transferred  to  the  waves.  The  model  has  been  extended  to  the  relativistic  and 
multi-mode  cases.  For  E = 0 the  nonrelativistic  results  are  in  good  agreement 
with  the  work  of  (OWM),  while  in  the  relativistic  case  the  model  reproduces  the 
nontrivial  features  of  the  computer  study  by  Lampe  and  Sprangle. 
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I.  INTRODUCTION 

One  of  the  more  striking  properties  of  hot  tenuous  plasmas  is  the  virtual 
disappearance  of  collisional  effects.  A well  known  consequence  of  this  behav- 
ior is  the  creation  of  runaway  electrons1  when  a sufficiently  strong  external 
DC  electric  field  is  applied  to  a plasma.  Although  various  aspects  of  this 
process  have  been  extensively  investigated,  there  has  been  a recent  resurgence2"4 
of  interest  in  this  topic.  The  reason  for  the  renewed  interest  can  be  attrib- 
uted to  the  successful  operation  of  toroidal  confinement  devices  capable  of 
entering  the  hot  plasma  regime  while  providing  sufficiently  long  confinement 
times  so  that  a runaway  population  can  be  established. 

The  study  of  the  properties  of  electron  runaways  is  of  interest  because 
it  provides  an  ideal  arena  for  testing  and  developing  various  theoretical  des- 
criptions of  fundamental  wave-particle  interactions  which  can  have  important 
consequences  beyond  the  narrower  realm  of  the  runaway  problem.  On  the  more 
practical  side,  a better  understanding  of  the  behavior  of  runaways  may  be  used 
as  a diagnostic  tool  for  uncovering  internal  phenomena  (e.g.,  magnetic  pertur- 
bations, transport,  instabilities)  not  directly  accessible  to  the  standard 
diagnostic  methods.  Also,  it  may  be  possible  that  in  certain  parameter 
regimes  the  large  streaming  energy  of  the  runaway  electrons  can  be  used  to 
generate  powerful  microwave  radiation,  or  even  to  achieve  the  direct  heating 
of  the  background  ions. 

The  motivation  for  the  present  study  has  been  provided  by  the  recent 
experimental  observations5  made  in  the  Microtor  tokamak  at  the  University  of 
California,  Los  Angeles.  By  monitoring  the  time  evolution  and  energy  of  the 
hard  x-rays  (>50  KeV)  emitted  from  this  device  one  deduces  the  maximum  energy 
of  the  runaway  population.  At  intermediate  values  of  the  plasma  density 
(ng  ~ 2 x 1012  - 1013  cm"3)  it  is  found  that  the  maximum  energy  of  the 
runaways  increases  monotonically  in  time  during  the  discharge,  eventually 
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reaching  levels  as  large  as  2-3  MeV.  However,  in  the  low  density  regime 
(ng  < 1012  cm'3)  it  is  found  that  the  maximum  runaway  energy  attained  holds 
constant  in  time  (5-10  msec)  during  a significant  fraction  of  the  discharge 
lifetime  (20-30  msec).  The  peak  energy  achieved  by  these  clamped  runaways 
is  of  the  order  of  300-400  KeV. 

Simultaneously  with  the  clamping  of  the  runaways  it  is  observed  that 
electromagnetic  radiation  in  the  neighborhood  of  the  ion  plasma  frequency, 

is  generated.  Spontaneous  emissions  near  the  electron  cyclotron  fre- 
quency, fle,  and  electron  plasma  frequency,  u>  , are  also  detected6  and  are 
present  in  both  density  regimes,  i.e.,  the  high  frequency  emission  is 
triggered  whether  or  not  the  runaways  are  clamped  in  energy.  However,  the 
lower  frequency  radiation  is  only  present  when  the  clamping  process  is  at 
work,  and  the  peak  of  this  emission  occurs  at  a frequency  which  scales 
proportionally  to  i/n^. 

Evidently,  such  an  interesting  behavior  in  the  complicated  environment 
of  a tokamak  plasma  can  stimulate  a wide  variety  of  explanations.  In  partic- 
ular, a simple  physical  picture  suggested  by  these  observations  and  whose 
interest  goes  far  beyond  the  previously  mentioned  experiment  is  the  following: 
a high  energy  electron  is  pushed  directly  by  the  external  DC  electric  field, 
however,  if  the  electron  can  resonate  with  a mode  supported  by  the  system 
then  a significant  part  of  the  push  may  be  transformed  to  the  form  of  wave 
momentum,  thus  resulting  in  the  clamping  of  the  particle.  This  type  of 
phenomenon  is  the  subject  of  the  present  study. 

Specifically,  this  study  is  concerned  with  the  evolution  of  a low  density 
beam  which  is  continuously  accelerated  by  an  external  DC  electric  field.  The 
beam  is  envisioned  to  meet  the  criteria  for  the  weak  cold  beam  description7* 8 
used  extensively  in  investigations  of  the  beam-plasma  interaction.  The  beam 
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is  assumed  to  propagate  along  a strongly  magnetized  plasma  capable  of  support- 
ing modes  whose  wavenumbers  are  quantized  by  the  geometry.  The  consideration 
of  such  modes  is  relevant  to  the  low  frequency  (u>  <<  w ) oscillations  in 
small  toroidal  machines  in  which  an  integral  number  of  wavelengths  fit  around 
the  torus.  A similar  situation  is  also  encountered  in  computer  simulations  in 
which  the  periodicity  length  introduces  the  quantization  of  the  Fourier  modes. 

The  mathematical  formulation  used  in  the  present  study  does  not  follow 
the  detailed  dynamics  of  the  beam  particles, which  would  require  a computer 
simulation.  The  method  used  here  relies  on  the  exact  conservation  equations 
for  the  total  energy  and  total  momentum  of  the  beam-wave  system.  To  bypass 
the  complexity  associated  with  the  simulation  of  beam  dynamics,  a WKB  model  is 
introduced.  The  foundations  of  this  model  can  be  traced  to  two  previous 
unconnected  investigations  whose  essence  is  briefly  described  in  the  following. 

Several  years  ago9  O'Neil  and  Malmberg  studied  the  topological  transition 
in  the  linear  dispersion  relation  for  the  cold  and  warm  beam-plasma  interaction. 
In  analyzing  the  transition  from  a cold  to  a warm  beam,  it  proved  useful  to 
consider  a two  parameter  Lorentzian  to  represent  the  beam  distribution  function. 
One  parameter  describes  the  average  velocity  of  the  beam,  while  the  other 
represents  its  thermal  spread.  The  results  obtained  with  this  parametrization 
of  the  problem  illustrated  clearly  the  continuous  transition  from  a hydrodynamic 
instability  to  a kinetic  Landau-type  of  interaction. 

The  other  investigation  which  has  a bearing  on  the  foundations  of  the  WKB 
model  used  here  concerns  a laboratory  experiment10  by  Starke  and  Malmberg.  This 
experiment  has  looked  in  detail  at  the  connection  between  the  well-known 
trapped  particle  sidebands  excited  by  a launched  large  amplitude  wave  and  the 
distortions  in  the  time  averaged  velocity  distribution  function.  The  remarkable 
outcome  of  this  experiment  is  that  the  behavior  of  the  spontaneously  excited 
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sidebands  can  be  explained  quantitatively  by  using  the  time  averaged  distribu- 
tion to  obtain  a local  linear  dispersion  relation. 

The  clues  that  have  been  extracted  from  the  two  previously  described 
results  to  solve  the  present  problem  are:  (1)  the  thermal  spreading  of  the 
runaway  beam  which  occurs  as  the  beam  encounters  a cavity  resonance  can  be 
reasonably  described  through  a two  parameter  Lorentzian;  the  two  parameters 
being  allowed  to  change  self-consistently  with  the  conservation  laws,  (2)  to 
bypass  the  calculation  of  the  dynamics  of  the  beam  particles,  one  uses  the 
slowly  changing  Lorentzian  into  the  linear  dispersion  relation  to  obtain  the 
time  evolution  of  the  wave  energy  and  momentum. 

To  check  the  validity  of  the  previously  described  formulation  we  have 
tested  its  predictions  against  the  detailed  particle  simulation  of  the  weak 
cold  beam  problem  by  O'Neil,  et  al.8  It  is  found  that  the  spatially  averaged 
technique11  introduced  in  this  work  predicts  the  correct  saturation  time;  the 
predicted  saturation  level  does  not  exhibit  the  trapped  particle  oscillations 
(as  expected)  but  rather  it  yields  an  average  value  which  ranges  between  the 
peaks  and  the  valleys  of  such  oscillations. 

When  the  effect  of  the  external  DC  field  is  retained,  one  finds  that  for 
certain  parameter  values  (e.g.,  the  amplitude  of  the  initial  noise  is  small) 
the  runaway  beam  can  sweep  through  a cavity  resonance  and  experiences  only  a 
slight  drag  and  a correspondingly  small  thermal  broadening.  However,  as  the 
level  of  the  initial  noise  is  increased  one  finds  that  the  beam  momentum  can 
be  clamped.  In  the  clamped  state  the  beam  momentum  increases  slowly  while 
the  amplitude  of  the  cavity  mode  grows  secularly  in  time.  In  this  regime  the 
external  push  by  the  DC  field  is  converted  into  wave  momentum  (DC  to  AC 
conversion) . 

L - - ■ mJ 
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The  spatially  averaged  formulation  has  been  extended  to  include  the  many 
mode  problem  as  well  as  the  relativistic  runaway  case.  The  clamping  behavior 
is  also  observed  when  these  effects  are  present.  The  predictions  of  the 
relativistic  formulation  with  the  DC  field  set  to  zero  have  been  compared 
with  the  computer  calculations12  of  Lampe  and  Sprangle.  It  is  found  that  one 
recovers  the  known  but  nontrivial  decrease  of  the  saturation  level  as  the 
parameter  s,  first  introduced13  by  Thode  and  Sudan,  is  increased.  In  addition, 
one  finds  that  the  actual  dependence  on  s is  in  excellent  agreement  with  the 
results  of  Lampe  and  Sprangle. 

The  investigation  of  the  restricted  runaway  beam  interaction  with  two 
modes  exhibits  a wealth  of  structure  that  uncovers  many  features  which  have 
been  previously  observed  in  the  laboratory10 » 1‘4> 1 5 in  the  study  of  trapped 
particle  sidebands.  Many  of  the  features  can  be  understood  in  terms  of  wave- 
particle-wave  interactions  in  which  the  beam  plays  the  role  of  an  intermediate 
agent. 

The  manuscript  is  organized  as  follows.  In  Sec.  II  the  problem  is 
defined  and  the  spatially  averaged  formulation  is  introduced.  Section  III 
discusses  the  application  of  the  formulation  to  the  weak  cold  beam  problem. 

Section  IV  presents  the  nonrelativistic  results  and  introduces  the  clamping 
behavior.  Section  V considers  the  relativistic  formulation.  Section  VI 
presents  a perturbation  theory  of  the  effects  produced  by  cavity  damping. 

Conclusions  are  presented  in  Sec.  VII. 
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II.  FORMULATION 

The  present  idealized  physical  model  consists  of  a low  density  beam 
(initially  cold)  which  is  pushed  by  an  external  DC  electric  field  Eq  that 
points  along  a strong  magnetic  field  Bq.  The  plasma  is  taken  to  be  a cylinder 
of  radius  a,  surrounded  by  a p,  rfectly  conducting  wall.  For  simplicity  the 
plasma  density  nQ  is  taken  to  be  radially  and  axially  uniform,  and  the  plasma 
cylinder  is  envisioned  to  be  axially  periodic,  with  periodicity  length  L. 

The  strong  guide  magnetic  field  points  along  the  axis  of  the  plasma  cylinder 
(z  axis). 

In  the  strong  magnetic  field  limit  (Bq  -*■“),  as  has  been  illustrated16 
by  Trivelpiece  and  Gould,  the  dispersion  relation  for  the  electromagnetic 
modes  (E  modes)  can  be  solved  analytically  and  is  obtained  from 

V2£z  + (k|  - k2)c.zEz  =0  (1) 

where  V2  refers  to  the  perpendicular  Laplacian  operator,  kQ  = u/c  corresponds 
to  the  vacuum  wavenumber  of  a signal  having  frequency  cu,  and  c is  the  speed 
of  light.  The  axial  component  of  the  dielectric  tensor  is  e . 

The  solution  of  Eq.(l)  takes  the  form 

Ez(z,r,0,t)  = AJn(k1r)exp{i (kz  + n6  - cot) } (2) 

in  which  r,0  represent  the  radial  and  angular  coordinates,  t is  the  time 
variable,  k is  the  axial  wavenumber,  A is  the  amplitude,  is  the  Bessel 
function  of  order  n,  and  n is  an  integer. 

The  finite  size  constraints  of  the  model  require  that 
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where  j is  an  integer  and  P refers  to  the  zeroes  of  J . Utilizing  these 
J ° nv  n b 

constraints  one  arrives  at  the  following  linear  dispersion  relation 

v 2 


(p  y 
nv 

"zz  ~ (kQa) 2-  (ka) 2 


(4) 


It  is  well  known  that  in  the  absence  of  the  beam  and  for  a cold  background 
plasma,  Eq.(4)  exhibits  two  branches:  (1)  the  fast  branch  consists  of  waves 


whose  phase  velocities  are  larger  than  c and  exists  for  frequencies  to  > o> 
(the  electron  plasma  frequency),  (2)  the  slow  branch  has  phase  velocities 


pe 


less  than  c and  exists  in  the  frequency  interval  to  < to  . It  is  this  latter 
branch  which  is  of  interest  to  the  present  problem  because  it  is  possible  for 
the  velocity  of  the  beam  to  coincide  with  the  phase  velocity  of  these  modes, 
thus  resulting  in  a strong  resonant  interaction. 

The  contribution  of  the  low  density  beam  to  the  finite  size  dispersion 
relation  enters  through  the  beam  susceptibility  X^,  where 


e = e°  + X, 
zz  zz  b 


(5) 


and  in  which  represents  the  contribution  of  the  background  plasma. 
Defining  the  generalized  finite  size  dielectric  of  the  system  e by 


e = e® 


(P  ) 
v nv' 


zz  (k^a )2-  (ka) 2 


(b) 


results  in  the  relation 


e(k,u>)  = - X 


(7) 


which  for  a low  density  beam  can  be  treated  in  analogy  with  the  small  cold 
beam  problem  , i.e.,  in  zeroth  order  E(k,w0)  = 0,  hence  the  perturbation 
produced  by  the  beam  can  be  described  through  to  = + Sin,  | Sin/uQ  | <<1 , 

and  is  determined  from 
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6u>  = - X. 


with  the  understanding  that  depends  functionally  on  6a>.  Once  the  depend- 
ence of  X^  on  6u  is  established,  Eq.(8)  can  be  used  to  obtain  5w,  whose 
imaginary  part  determines  the  growth  or  damping  of  the  cavity  mode  supported 
by  the  background  plasma.  In  the  present  work  Eq.(8)  is  used  in  its  general- 
ized WKB  sense,  i.e.,  one  calculates  the  slow  time  evolution  of  X^  = X^(t,6w). 

In  general  to  find  X^  one  needs  to  follow  the  detailed  dynamics  of  the 
beam  electrons.  To  extend  such  an  exercise  beyond  the  linear  regime  one 
typically  resorts  to  a digital  computer  in  order  to  calculate  the  complicated 
trapped  particle  orbits.  However,  in  the  present  study  we  wish  to  avoid  such 
a procedure  by  relying  on  a spatially  averaged  description  based  on  the  exact 
conservation  laws  for  the  beam-wave  system. 


The  derivation  and  meaning  of  the  conservation  laws  are  easier  to  under- 
stand if  one  represents  the  beam  in  terms  of  N discrete  particles.  Each 
particle  is  labeled  by  the  index  l,  1 < l < N,  and  in  the  nonrelativistic 
limit  their  motion  is  determined  by 


m dt 


eEQ  - e { E(t)Jn(k1rJl)exp[i(kzJl  + ne^-u^t)]*  c.c.}  (9) 


in  which  e,  m refer  to  the  electron  charge  and  mass,  v is  the  axial  component 
of  the  velocity  of  particle  l with  coordinates  (z^,  0^,  r ) , and 


t 

!(t)  = E (O)exp  - i dt'6uj(t') 
z 

' o J 


with  6w  complex. 


After  Fourier  analyzing  in  time,  Maxwell's  equations  take  the  form 


v x v x r - i ■ + (— ) 

* C ^ 
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where  j refers  to  the  oscillating  plasma  current,  which  in  the  strong  Bq  limit 
points  essentially  in  the  z direction.  Accordingly, 

2 

(V  x v x E)z  = - V1Ez  ♦ 3z(V1  • EJ  (12) 

and 


3zEz  + V1  • Ex  = 4irp 


where  p refers  to  the  oscillating  electron  charge  density, 
and  (13)  in  Eq.(ll)  yields 


i2p  _i_  c _ a — ^ _ i4TTU) 


V.E  + 3ZE  + (— ) E = 4tt3  p 

1 7 Z Z 'p/  7 7. 


Operating  on  both  sides  of  Eq.  (14)  with 


(13) 

Using  Eqs. (12) 


(14) 


L/2  r 2n 

l - J ^7  exp[-i(kz  + n6) ] 

-L/2 


(15) 


and  using  the  conservation  of  charge  constraint  together  with  the  periodicity 
requirement  along  z results  in 


VfE  + (kz  - kz)E  = 
1 z o z 


4u 


- 1 r <ko  • k2)» 


(16) 


where  the  bar  notation  over  Ez>  p refers  to  the  Fourier  transformed 
quantities  over  z and  6. 

The  quantity  p contains  the  contribution  pq  arising  from  the  background 
particles  (which  are  assumed  to  behave  linearly),  as  well  as  the  contribution 
p"b  due  to  the  beam  particles,  i.e., 


. . (0  ? 
F-  £ <-K)  2e 

4tt  ' a)  2 


♦ Pv 


(17) 
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Inserting  Eq.(17)  into  Eq.(16)  and  integrating  over  the  radial  coordinate 
with  the  appropriate  weight  functions  results  in 

r3  a V,¥ 

2 ^ J (k.r)  { tt5 — ♦ c°  E } 

a2  nv  * ’ k2  - kz  zz  z 


IN 

_ 4ne  ,V  l J [k1r,(t)]exp{-i[kz  (t)  ♦ n9  (t)]} 

1 “IT  '•N  J t=l 

where  n^  represents  the  density  of  beam  electrons. 

For  low  density  beams  (i.e.,  n^/n^  <<:  1)  the  left-hand  side  of  Eq.(18) 
vanishes  to  zero-order  in  Sid,  thus  yielding  the  solutions  described  by  Eqs.(2) 
and  (4).  To  account  for  the  interaction  between  the  beam  and  the  cavity  mode 
one  must  retain  the  first  order  contributions  in  Sid  which  arise  in  the  left- 


hand  side  of  Eq . (18) , namely 


a,  I (2k  /c)  9e° 

— t J (k.r)  - ■-■■■  y-.-A-'-A  VfE  + — — F. 

a2  nv  t (k2-k2)2  l z 3<d  z 


where  it  is  understood  that  the  terms  inside  the  brackets  are  evaluated  at 
id  = <dq,  i.e.,  the  unperturbed  frequency  of  the  cavity  mode.  Since  the  zero 
order  Eqs.(l)  and  (4)  are  satisfied  identically,  the  expression  in  Eq. (19) 


takes  the  form 


2(k  a)2  /e°  I2  h 


zz  zz 


P^T  [J2+1(kia)]SU(t)E(t) 


Using  the  generalized  dielectric  c defined  in  Eq.((>)  one  finally  arrives  at 
the  equation  governing  the  time  evolution  of  E(t) 


(-^0  fJ2  , (k , a)  ] 5-E  = — ~ ) 
v3id'  1 n+lv  1 'J  dt  Nk  (r 

o t 


l Jn(kir^)exp{-i(kzjl+nejl-iDot)}  (21) 


To  obtain  the  time  evolution  of  the  total  momentum  of  the  system  one 


returns  to  Eq. (9)  and  sums  over  l,  making  use  of  Eq.(21),  to  yield 
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at  {"4  (■*!£}  - «V„  (22) 

which  states  that  the  total  momentum  of  the  beam-wave  system  is  increased  by 
the  presence  of  the  external  DC  electric  field. 

In  calculating  the  conservation  of  energy  relationship  it  proves  useful 
to  introduce  the  relative  coordinate  C0  = z.  - (w  /k)t  and  the  associated 
velocity  l . Multiplying  both  sides  of  Eq.(9)  by  and  summing  over  all 
beam  particles  yields 


d v 

At  1 7 


eEjk,  - eKE(t)5j;Jn(kire)exp[-i(k5jl+n0a)]fc.c 
1 1 1 


but  in  the  strongly  magnetized  limit 


IiJlexp[i(kCil+n0ll)]Jn(k1ri,)  = £ Jn(kjlr)l)exp[i  (k^+n©^  ] ' (24) 


using  Eq .(21)  this  expression  becomes 


-L_  5L_  r Nk_.  j2  ^ g*i 

ik  dt  Mnen,  ^3oj'  ' dt 

b 

Inserting  Eq.(25)  in  Eq.(23)  and  recognizing  that 

E 3?2E*  - E*  E = - Ht  [2ReC6u>)  |E|2] 


where  Re  refers  to  the  real  part,  results  in 


it  %l  t * 2ReW“)(lz)Jn.i(kia'[Ji!r, 


)k  N o 


which  states  that  the  time  rate  of  change  of  the  total  energy  of  the  beam- 
wave  system  is  equal  to  the  power  delivered  by  the  DC  electric  field. 


M 
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Equations  (22)  and  (23)  are  the  exact  conservation  laws  for  the  runaway 
beam-cavity  interaction  used  in  the  present  work.  Similar  expressions  have 
been  found  previously  in  the  discussion  of  the  weak  cold  beam  problem,  both 
in  the  nonrelativistic8  and  relativistic17  cases,  and  also  in  the  calculation 
of  the  nonlinear  frequency  shift18  due  to  particles  trapped  by  an  electron 
plasma  wave.  In  the  latter  study  it  has  been  shown  that  to  lowest  order,  the 
conservation  of  momentum  regulates  the  damping  or  growth  of  the  wave,  while 
the  conservation  of  energy  requirement  produces  phase  changes  in  the  wave.  In 
the  present  work  both  of  these  effects  play  a crucial  role  in  determining  the 
evolution  of  the  runaway  beam.  For  completeness,  it  should  be  mentioned  that 
quasilinear  theory  retains  only  the  effects  associated  with  momentum  changes 
and  completely  neglects  the  phase  changes  introduced  by  Eq.(27). 

As  can  be  seen  from  Eqs. (22)  and  (27)  the  conservation  relations  depend 
on  two  spatially  averaged  properties  of  the  beam,  i.e.,  the  average  velocity 

/N,  and  the  average  of  the  squares  of  the  velocity  £(£.)2/N.  If  one  finds 

l * 

the  time  evolution  of  these  quantities  then  the  behavior  of  | E | 2 and  6u>  can 
be  determined.  Such  a procedure  can  be  implemented  if  the  orbits  of  the  beam 
particles  are  known.  In  the  case  of  a launched  large  amplitude  wave  this 
method  of  calculation  has  been  implemented  analytically,18’18  while  in  the 
small  cold  beam  problem  the  orbits  have  been  calculated  numerically. 8 In  the 
present  work  we  are  interested  in  using  the  conservation  equations  in  the 
reverse  order.  Essentially  the  idea  is  to  determine  6w  and  | E | 2 first,  and 
then  use  the  conservation  laws  to  extract  the  previously  mentioned  averages. 

Of  course,  such  a procedure  must  be  self-consistent  thus  implying  that  a 
method  of  closure  must  be  implemented. 

In  order  to  calculate  6<u(t)  from  Eq.(8)  one  needs  to  know  X^(t,6w). 
Physically,  one  expects  that  as  the  runaway  beam  passes  through  a cavity 
resonance  its  average  velocity  is  altered  and  simultaneously  it  acquires  a 
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certain  thermal  spread.  This  implies  that  the  initial  cold  beam  is  slowly 
transformed  into  a warm  beam  which  suffers  a recoil  due  to  the  fact  that  it 
excites  a cavity  mode.  Thus  one  must  deal  with  the  transition  from  the  hydro- 
dynamic  interaction  (two-stream)  to  the  kinetic  interaction  (Landau  type). 
O'Neil  and  Malmberg  have  investigated  the  details  of  such  a transition  in 
the  purely  linear  regime,  and  have  found  that  it  is  extremely  useful  to 
parametrize  the  transition  by  means  of  an  equivalent  beam  succeptibility 


J* 


(us-kvb+ikv)2 


(28) 


where  w2b  = wpe(nt/n0)>  v^  is  the  average  velocity  of  the  beam  and  v repre- 
sents the  average  thermal  spread.  In  the  limit  v = 0,  Eq.(29)  reduces  to 
the  cold  beam  result,  while  for  v / 0 this  succeptibility  can  be  identified 
with  a beam  having  a Lorentzian  distribution. 

The  closure  procedure  consists  of  using  Eq.(28)  in  its  WKB  sense,  i.e., 
v^=  vb(t)  and  v = v(t),  thus  describing  a slowly  changing  beam  which  under- 
goes a transition  from  cold  to  warm  as  it  interacts  with  the  cavity  mode. 

To  make  the  scheme  self-consistent  one  relates  vb  and  v to  the  spatially 

averaged  quantities  appearing  in  the  conservation  laws,  namely 

1/2 

_ . : _ . ' 

(29) 


; v(t)  = i 

r 2 

l—  - 1 
h n 

1 y ^ ' 
o N 

2 

l 

l4  ) 

With  Eq . (28)  the  WKB  (in  time)  dispersion  relation  takes  the  form 

. w2  (n,  /n  ) 

, 3c^P  pe  b o 

t3o)}l5w(t)  " [6u>(t)-kvb(t)  + ik7(t)]2 


(30) 


and  with  the  understanding  that  vb  is  measured  relative  to  the  wave  frame. 
Introducing  the  following  scaling 
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W = 6(J)/S1,  X = flt 


p = kvb/«,  S = kv/J5 


u = 


k2(fS)Jn+l(kia)lEl2  /C4irmnbfi) 


= LIQ 


kCn'  2 


2N  i « 


£=  (ekEo)/(mn2) 


one  obtains  the  scaled  equations  of  the  problem 


(31) 


37  (P  + u)  = € 

(32) 

[K  + 2Re(w)u]  = <fp 

(33) 

w(w  - p ♦ iS) 2 = 1 

(34) 

S = \ (2K  - p2)*5 

(35) 

The  simplification  achieved  by  the  spatially  averaged  description  becomes 
evident  by  integrating  Eqs.(32)  and  (33)  directly 


f — \ 

H 

II 

H 

1 

f — \ 

H 

1 

[p(0)  ♦ u(0) ] 

(36) 

K(x)  = 2Re[w(x)]  ♦ 

fdx'p(x')  - 2Re[w(0)]u(0) 

J 0 

(37) 

with 


U(T) 


u(0)exp{2 


ft 


dx' Im[w(T')] ) 
0 


(38) 


and  realizing  that  w(x)  can  also  be  obtained  in  closed  analytic  form  from 
Eq . (34)  by  means  of  the  well-known70  formula  for  the  roots  of  a cubic 
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equation  (not  worth  writing  out  in  detail  here).  The  latter  simplification 
follows  from  the  choice  of  Eq.(28),  other  methods  of  closure  yield  more  compli- 
cated expressions.  In  practice  the  solution  to  the  present  problem  reduces 
essentially  to  solving  a set  of  coupled  algebraic  equations  which  can  easily  be 
done  with  very  little  computing  effort. 

It  should  be  noted  that  the  procedure  used  in  generating  Eqs. (32)-(35)  can 
be  implemented  for  many  modes.  In  this  more  general  case  the  resulting  equations 
are  slightly  changed  to  the  form 


where  (3e/3u))^  refers  to  (3e/3w)  evaluated  at  the  unperturbed  frequency  of 
mode  k,  i.e.,  It  is  clear  that  for  k = the  expressions  in  Eq.(43) 

reduce  to  those  of  Eq.(31). 


r 
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III.  THE  WEAK  COLD  BEAM  SATURATION 

It  is  clear  that  the  spatially  averaged  description  neglects  a variety 
of  single  particle  dynamical  effects  which  have  been  found  to  play  an  impor- 
tant role  in  previous  investigations7* 8 of  the  beam-plasma  interaction. 
Outstanding  among  these  is  the  trapping  of  the  beam  by  the  fastest  growing 
linearly  unstable  mode.  The  trapping  effect  saturates  the  hydrodynamic 
instability  and  causes  periodic  oscillations  of  the  wave  amplitude.  Physic- 
ally, the  trapping  process  causes  the  cold  beam  to  acquire  a velocity  spread 
which  in  turn  stops  the  growth  and  produces  the  temporary  damping  of  the 
wave.  However,  due  to  the  strong  coherency  of  this  system  the  velocity 
spread  is  only  partially  randomized.  Accordingly,  a significant  fraction 
of  the  fake  heating  is  reversible  and  manifests  itself  in  the  regrowth  of 


the  wave.  This  sequence  of  events  repeats  cyclically  until  the  system 
phase  mixes  or  other  processes  occur  which  make  the  beam  heating  practically 
irreversible. 

The  spatially  averaged  description  can  handle  the  effects  associated 
with  the  increase  in  the  velocity  spread.  However,  the  model  intrinsically 
assumes  that  the  velocity  spread  is  an  irreversible  process.  Consequently, 
the  model  can  not  reproduce  the  trapped  particle  amplitude  oscillations, 
instead  it  gives  an  answer  which  is  an  average  over  these  oscillations. 

To  check  the  predictive  power  of  the  model  we  consider  the  limit  in 


which  there  is  no  external  DC  electric  field  present,  and  at  first  consider 
a single  mode.  In  this  limit  the  problem  reduces  to  the  weak  cold  beam- 
plasma  interaction  investigated  by  Shapiro,  et  al.7and  by  O'Neil,  et  al.8 

Figure  1 shows  the  comparison  between  the  spatially  averaged  model 
(continuous  curve)  and  the  results  of  O'Neil,  et  al.8  (dashed  curve). 


In  the  upper  portion  of  Fig.  1 one  plots  the  time  evolution  of  the  scaled 
wave  amplitude  */u.  It  is  seen  that  the  two  results  overlap  in  the  linear 


regime  and  reach  saturation  at  the  same  time.  As  expected,  after  saturation 
the  spatially  averaged  model  yields  a constant  amplitude,  while  the  computer 
simulation  predicts  amplitude  oscillations  about  the  saturation  level  of  the 
averaged  model.  On  the  bottom  of  Fig.  1 one  observes  the  time  evolution  of 
the  scaled  average  momentum  of  the  beam,  and  shows  that  the  beam  slows  down 
as  the  instability  saturates. 

Although  the  beam  momentum  is  converted  into  wave  momentum,  as  illus- 
trated in  Fig.  1,  this  effect  is  not  the  primary  reason  for  the  saturation. 

To  isolate  the  physics  behind  the  saturation  we  have  considered  the  evolution 
of  the  system  when  one  legislates  that  S = 0 for  all  t,  i.e.,  one  does  not 
allow  the  beam  to  acquire  a thermal  spread  as  it  recoils.  The  resulting 
behavior  is  shown  in  Fig.  2,  where  one  observes  that  the  instability  does 
not  saturate  in  this  case.  This  result  shows  that  it  is  crucial  to  retain 
the  thermal  spreading  of  the  beam  in  order  to  correctly  describe  the  satura- 
tion illustrated  in  Fig.  1. 

The  proper  physical  interpretation  of  the  saturation  is  that  as  the 
beam  acquires  a thermal  spread  the  hydrodynamic  instability  evolves  into 
the  kinetic  instability.  The  consequence  of  this  transition  is  illustrated 
in  Fig.  3,  where  one  observes  the  dependence  of  the  growth  rate  r on  p for 
a cold  beam,  S = 0,  and  for  a slightly  warm  beam,  S = 0.5.  It  is  seen  from 
Fig.  3 that  as  the  beam  acquires  a thermal  spread  the  region  of  unstable 
behavior  shrinks  and  moves  toward  larger  values  of  p.  However,  the  conser- 
vation of  momentum  constraint  requires  that  the  momentum  p must  decrease  as 
a function  of  time.  Thus,  as  the  beam  recoils  it  eventually  passes  through 
the  marginally  stable  point  (i.e.,  r = 0)  bringing  about  the  saturation. 

For  completeness  we  exhibit  in  Fig.  4 the  self-consistent  evolution 
of  the  kinetic  energy  K,  thermal  spread  S,  and  the  frequency  shift  ft  = Re(w) 


A 


corresponding  to  the  solid  curve  of  Fig.  1.  It  should  be  noted  that  the 
kinetic  energy  is  measured  relative  to  the  initial  wave  frame,  hence  the  beam 
recoil  shows  up  as  an  increase  in  K. 

Another  topic  which  can  be  investigated  by  the  formalism  is  the  genera- 
tion of  sidebands  due  to  wave-particle-wave  interactions.  This  interaction 
has  been  recently  identified  as  the  fundamental  process  governing  the  excita- 
tion of  the  well-known  trapped  particle  sidebands.  The  essence  of  this 
interaction  is  that  a wave  having  a certain  phase  velocity  can  cause  a non- 
local distortion  in  the  spatially  averaged  velocity  distribution  function, 
thus  modifying  the  damping  or  growth  of  other  modes  having  different  phase 
velocities.  Such  a process  can  be  present  in  the  beam-plasma  interaction,  as 
is  illustrated  in  Fig.  5. 

In  Fig.  5 we  exhibit  the  time  evolution  of  the  amplitude  of  two  modes, 

^ obtained  from  the  self-consistent  numerical  solution  of  Eqs . (39) - (42) 

for  a fixed  choice  of  k0/k ^ = 0.7,  and  for  different  values  of  dispersion, 
i.e.,  6,  = 2.0,  1.5,  0.5.  The  scaling  wavenumber  k^  corresponds  in  this 

case  to  the  fastest  growing  mode  in  the  linear  stage  of  the  beam-plasma 
instability.  It  is  seen  in  Fig.  5 that  early  in  time  both  modes  exhibit 
growth  for  all  values  of  6.,.  This  is  the  linear  stage  during  which  the  beam 
is  not  significantly  perturbed.  However,  in  the  nonlinear  regime  one 
observes  different  responses  for  different  values  of  6.,.  For  6_,  = 2.0  one 
has  the  more  familiar  situation  arising  due  to  the  heating  of  the  beam  by 
the  growing  waves.  As  the  beam  acquires  a thermal  spread  (mainly  due  to  the 
growth  of  k^)  mode  k^  undergoes  the  transition  from  hydrodynamic  to  kinetic 
behavior,  thus  resulting  in  its  damping.  The  final  state  consists  of  the 
saturated  fastest  growing  mode  kj. 
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It  is  seen  in  Fig.  5 that  as  the  dispersion  is  decreased  to  = 0.5, 
the  fastest  growing  mode  undergoes  a continuous  transition  that  results  in 


its  eventual  decay,  i.e.,  a sideband  decay  whereby  the  energy  in  mode  k 


1 


is  transferred  to  the  neighboring  mode  k^.  This  is  the  type  of  decay 
responsible  for  the  excitation  of  the  trapped  particle  sidebands1 0 _1 h and 
should  be  noted  that  it  has  nothing  to  do  with  the  parametric  couplings 
previously  invoked21  to  explain  such  a behavior. 

To  better  understand  the  reason  behind  the  nonlinear  decay  of  the 

fastest  growing  mode  for  = 0.5,  one  should  consider  the  self-consistent 

modification  in  both  the  real  and  imaginary  parts  of  Wj  and  w^.  The 

relevant  time  evolution  of  n.  = Re(w.)  and  r . = Im(w.)  is  shown  in  Fig.  6 

3 3 3 3 6 

for  6^  = 0.5.  Early  in  time  is  slightly  larger  than  and  j S2 ^ | < | £2^ | . 

As  the  beam  acquires  a thermal  spread,  decreases  to  zero  as  it  would 
normally  do  just  before  saturation  for  the  single  mode  problem.  However, 
because  mode  k^  is  now  present  the  beam  continues  to  recoil  and  spread, 
thus  causing  mode  kj  to  overshoot  the  marginally  stable  point  (Tj  = 0)  and 
eventually  enter  the  region  of  damping.  As  mode  kj  damps,  the  recoil  of 
the  beam  is  slowed  down  thus  preventing  mode  k^  from  crossing  its  marginally 
stable  point  (i.e.,  F^  = 0) . Also,  from  the  scaling  of  Eq .(41)  one  observes 
that  for  a fixed  thermal  spread  S,  mode  k_,  samples  an  effective  thermaliza- 
tion  which  is  0.5  smaller  than  that  seen  by  k..  As  a consequence, 
suffers  a much  greater  change  from  its  linear  value  than  as  is  seen  on  the 
bottom  of  Fig.  6. 

The  present  model  can  describe  also  the  cascading  of  the  spectrum  to 


lower  phase  velocities  (higher  k) . This  effect  is  expected  to  be  responsible 
for  the  long  term  slowing  down  of  the  beam  and  its  final  merger  with  the 
distribution  function  of  the  background  plasma.  To  isolate  this  effect  we 
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consider  two  modes  of  widely  spaced  phase  velocities  = 3.0,  k^/k^  = 0.7, 
and  legislate  that  the  initial  beam  velocity  lies  between  them,  i.e., 
p(0)  = 2.0.  In  this  case  k^  corresponds  to  a mode  which  is  initially  stable 
while  k0  is  slightly  unstable,  as  predicted  from  Fig.  3 for  S = 0.  The 
corresponding  evolution  of  the  wave  amplitudes  is  illustrated  on  the  top 
part  of  Fig.  7,  where  it  is  seen  that  k^  grows  initially,  saturates  early  in 
time  and  eventually  decays.  As  mode  k.,  decays  one  observes  that  mode  kj  is 
excited.  The  reason  for  the  late  excitation  of  mode  k^  is  that  the  satura- 
tion of  k,  causes  the  beam  to  recoil,  as  indicated  on  the  bottom  portion  of 
Fig.  7.  As  p is  reduced,  mode  k^  enters  the  region  of  growth  (as  in  Fig.  3), 
thus  bringing  about  the  cascading  to  lower  wavenumbers. 
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IV.  CLAMPING  BEHAVIOR 

We  proceed  to  consider  the  behavior  of  the  beam  when  an  external 
DC  electric  field  is  present.  In  addition  to  the  strength  of  the  DC  field 
one  needs  to  legislate  at  t = 0 the  momentum  of  the  beam  p(0),  its  initial 
velocity  spread  S(0),  and  the  initial  amplitude  of  the  various  modes,  u^(0). 

Figure  8 illustrates  the  time  evolution  of  the  wave  amplitude  and  beam 
momentum  for  a case  in  which  a cold  beam  interacts  with  a single  mode  having 
a relatively  small  initial  amplitude,  i.e.,  u^ (0)  = 1.0  * 10-3,  p(0)  = - 5.0, 
S(0)  = 0,  £*  = 1.0.  It  is  seen  in  Fig.  8 that  for  these  parameters  the  beam 
is  able  to  run  through  the  phase  velocity  of  the  cavity  resonance  while 
experiencing  only  a slight  recoil.  This  small  recoil  arises  because  for 
p < 0 the  beam  drives  the  cavity  mode  unstable  (as  in  Fig.  3).  However, 
since  the  beam  is  continuously  accelerated  by  the  external  DC  field  its 
velocity  eventually  becomes  larger  than  the  phase  velocity  of  the  cavity 
mode,  hence  shutting-off  the  instability.  For  small  initial  noise  levels 
and  large  DC  fields  the  process  consists  of  a runaway  beam  radiating  energy 
into  a cavity  mode  for  a finite  time  and  proceeding  to  increase  its  velocity 
until  it  encounters  a new  resonance  with  another  mode,  thus  repeating  the 
cycle.  This  behavior  would  show  up  experimentally  in  the  form  of  a cascade 
in  the  spectrum  of  the  observed  radiation  toward  lower  frequencies. 

An  interesting  phenomenon  occurs  when  either is  decreased  or  u^(0)  is 
increased.  By  decreasing  ^ one  increases  the  time  of  interaction  between 
the  runaway  beam  and  the  cavity  mode,  while  by  increasing  u^fO)  the  absolute 
amplitude  of  the  wave  can  be  increased.  Both  of  these  changes  can  make  the 
rate  of  momentum  transfer  to  the  wave  comparable  to  the  magnitude  of  £ , thus 
preventing  the  beam  from  running  away  in  velocity.  Such  a clamping  behavior 
is  illustrated  in  Fig.  9 for  p(0)  = - 5.0,  S(0)  = 0,  £"  = 1.0,  as  in  the  case 
of  Fig.  8 but  now  Uj(0)  = 5.0  * 10~3.  It  is  seen  in  Fig.  9 that  early  in 


time  the  beam  momentum  exhibits  the  typical  free-fall  behavior  associated 
with  runaways.  However,  for  p a 0 the  beam  recoils  and  its  momentum  proceeds 
to  grow  at  a much  slower  rate  while  the  wave  amplitude  (and  wave  momentum) 
continue  to  increase  secularly.  The  clamping  behavior  exhibited  in  Fig.  9 
is  an  example  whereby  an  external  push  on  a material  particle  can  be  trans- 
formed into  field  momentum.  Alternatively,  one  can  say  that  the  beam  acts 
as  a converter  of  PC  to  AC  energy. 

The  top  of  Fig.  10  illustrates  the  time  evolution  of  the  effective  force 
acting  on  the  beam.  Early  in  time  this  force  is  equal  to  the  DC  push  5*  = 1.0, 
and  later  on  it  is  seen  that  the  effective  force  decreases,  reverses  sign  for 
a short  time,  and  finally  it  settles  at  the  reduced  value  of  5 * 10-3.  One 
sees  in  the  bottom  of  Fig.  10  the  evolution  of  the  kinetic  energy  K and  the 
thermal  spread  S during  clamping. 

The  functional  dependence  of  the  slow  time  evolution  observed  after  the 
clamping  stage  sets  in,  illustrated  in  Figs.  9 and  10,  can  be  seen  to  be  a 
direct  consequence  of  Eq. (39)  in  the  regime  of  small  growth  rate  6,  but  large 
amplitude  so  that  u(x)  * u(l  + 20t)  with  u defined  by  the  condition  that 
dp/dx  = 0,  i.e.,  2Tu  The  slowly  changing  momentum  and  kinetic  energy 

are  of  the  form 


P - P + 2(6  -r  )ut 

K = K ♦ (£p  + 2wu)  x t2 


(44) 


where  the  bar  notation  denotes  constants  evaluated  just  after  clamping  is 
attained,  and  0 - T is  a small  quantity. 

In  estimating  the  threshold  value  of  the  Uj (0)  required  to  obtain 
clamping  for  a given  p(0)  and  C it  is  useful  to  consider  the  prediction 


made  by  the  first  Born  approximation 


ui(°)  >£exp  - (rM  ♦ r0) 


in  which  r refers  to  the  maximum  linear  growth  rate  and  r to  the  growth  rate 

M O 


for  p (0)  (as  predicted  from  Fig.  3).  By  examining  the  clamping  behavior  for 


a wide  range  of  initial  conditions  we  have  found  that  Eq.(45)  gives  good 


results  for  intermediate  values  of  p(0),  i.e.,  p(0)  = - 5.0.  As  an  illustra- 


tion, for  p(0)  = - 5.0,  1.0,  Eq.(45)  predicts  Uj(0)  > 1.6  x 10~3  and  by 

searching  in  parameter  space  one  finds  that  1.0  * 10"3  < u^  (0)  < 4.0  * 10'3. 


However,  Eq.  (45)  fails  to  predict  the  correct  threshold  for  higher  and  lower 


values  of  p,  i.e.,  p(0)  < - 7.5  and  - 2.5  < p(0).  The  physics  behind  the 


failure  is  directly  related  to  the  thermal  spreading  of  the  beam  and  plays 


opposite  roles  for  large  and  small  values  of  |p(0) 


For  large  | p (0) | the  beam  acquires  a finite  velocity  spread  as  it 


approaches  the  resonance,  thus  the  effective  growth  rate  is  reduced,  hence 


requiring  a larger  value  to  achieve  clamping.  For  small  |p(0)j  the  beam 


penetrates  deeper  into  the  p > 0 region.  However,  for  this  situation  the 


thermal  spread  elongates  the  region  of  growth  toward  larger  p (as  in  Fig.  3) 


thus  implying  that  the  effective  interaction  time  with  the  cavity  mode  is 


longer  than  estimated  by  the  first  Born  approximation.  As  a consequence,  the 


threshold  Uj(0)  is  lower  than  predicted  by  such  a method.  An  illustration  of 


the  deeper  penetration  at  smaller  p(0)  is  illustrated  in  Fig.  11  for 


p(0)  = - 2.5,  S(0)  = 0,  = 1.0,  Uj (0)  = 5.0  * 10'3.  It  is  seen  in  Fig.  11 


that  the  beam  penetrates  up  to  p = 0.8  (to  be  compared  with  Fig.  9)  and  then 


recoils  to  the  region  p < 0 while  the  wave  amplitude  grows  secularly. 


The  clamping  process  can  also  take  place  when  more  than  one  mode  is 


considered.  Figure  12  exhibits  the  behavior  for  two  neighboring  modes  with 


k2/kj  = 0.7,  6 2 = 3.0,  u2(0)  = Uj(0)  = 1.0  * 10'2,  p(0)  = - 5.0  and  = 1.5. 


1 


. . 


On  the  top  part  of  Fig.  12  one  sees  the  time  evolution  of  the  wave  amplitudes, 
and  it  is  observed  that  the  slower  mode  saturates  earlier  and  attains  a 
lower  amplitude  level  than  mode  k9.  For  large  t it  is  observed  that  both 
modes  exhibit  the  characteristic  secular  growth  associated  with  the  clamping 
of  the  beam.  By  examining  the  evolution  of  the  beam  momentum  shown  on  the 
bottom  of  Fig.  12,  one  finds  that  the  beam  actually  runs  through  the  phase 
velocity  of  mode  kj  until  it  encounters  mode  k.,,  which  clamps  the  beam  as  in 
the  single  mode  example  of  Fig.  9.  By  setting  u^fO)  = 0 one  verifies  that 

mode  k2  is  responsible  for  this  effect,  and  in  such  a situation  one  finds 

that  mode  k^  does  not  exhibit  the  late  time  secular  growth;  it  simply  saturates. 
The  dominant  multi-mode  effect  on  the  clamping  process  is  that  in  addition  to 
exciting  the  main  mode  responsible  for  the  slowing  down  of  the  beam,  neigh- 
boring modes  are  driven  unstable  due  to  the  Landau-type  of  growth  associated 
with  a beam  having  a finite  velocity  spread. 

The  effect  of  finite  velocity  spread  on  the  multi-mode  interaction  is 
illustrated  in  Fig.  13  for  S(0)  = 1.0,  u.,(0)  = Uj(0)  = 5.0  * 10-2,  with  all 

other  parameters  as  in  Fig.  12.  It  is  seen  on  the  top  part  of  Fig.  13  that 

both  modes  exhibit  an  early  damping  phase  followed  by  growth  and  eventual 
saturation.  This  behavior  is  an  entirely  kinetic  effect  associated  with  the 
sampling  of  the  negative  slope  of  the  velocity  distribution  function  (i.e., 

3f/3v  < 0)  before  the  beam  reaches  the  resonance.  Following  this  stage  there 
is  a period  of  growth  when  the  3f/3v  > 0 part  of  the  beam  passes  through  the 
cavity  modes  in  velocity  space,  as  suggested  on  the  bottom  part  of  Fig.  13. 
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V.  THE  RELATIVISTIC  PKOBLEM 

The  relativistic  runaway  beam  problem  can  also  be  handled  through  the 
appropriate  conservation  laws  of  the  beam-wave  system.  These  laws  can  be 
derived  in  analogy  with  the  nonrelativistic  case  and  for  the  single  mode 
interaction  take  the  form 


d 

dt 

d_ 

dt 


]_21 

w 


= en  E 

b o 


(46) 


mc2y  , | F | 2 

-V  - \lpi  * 


- vQ) 


= en,  E T 

b o ■*  N 


(47) 


where  P is  the  momentum  of  the  l beam  particle  and  y is  the  corresponding 


relativistic  factor.  In  Eq.  (47)  v^  is  the  phase  velocity  of  the  cavity 
mode  and  vfc  = (P^/m) [ 1+ (P^/mc) 2] " s. 

The  dispersion  relation  for  this  problem  is 


6W  = 


pb 

y3  [<$w-k(v-vo)  + ikA/  (my^)  ] 2 


(48) 


where  v refers  to  the  average  velocity  of  the  beam  and  A represents  a small 

momentum  spread,  i.e.,  A/mv  <<  1.  The  full  frequency  of  the  mode  is 

u = kv  +6w. 
o 

The  exact  conservation  laws  shown  in  F.qs.(46)  and  (47)  can  be  consid- 
erably simplified  if  one  chooses  to  describe  the  evolution  of  the  beam  sfter 
it  has  acquired  a large  momentum  Pq  relative  to  the  lab,  so  that 

yl  = Yo  (1  + “V*  p = Po(1  + aY0qP  (49) 


with  the  small  parameter  being 


<o2. 

J*L 

(— ) 

LWJ 


1/3 


a 


(50) 
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and  yo  is  obtained  from  Pq  through  the  relation  y(P)  = [1+ (P/mc) 2]*s.  This 
relation  also  provides  a connection  between  n0  and  q for  the  limit  a « 1, 

Xj  K 


namely 


|Vo|2,1(  l(Pol2  2 

\ ' — 72  * 2 K “>e 

' o 1 


which  implies  that 


k(vi  - V 


Y0°“U 


Introducing  the  natural  scaling  of  the  problem  w = 6w/a,  S = kA/ (may3) , 
transforms  Eq.  (48)  into 


r . iyo 

W W + 1 

L W 


S - q = [1  + an]'3 


with  q - Iq^/N,  n = £nf/N.  Working  to  lowest  order  in  q results  in  the 
dispersion  relation 


w [ w + iS Cl  + sq)'3  - q]2  = (1  + sq)'3 

which  is  seen  to  depend  explicitly  on  the  value  of  the  parameter 
s = (a/kvQ) (vo/c) 2Y0"2  which  was  originally  introduced13  by  Thode  and 
Sudan  to  study  the  beam-plasma  interaction. 

Using  Eq.(49)  in  Eq.(47)  and  scaling  t = at  yields 

dU  _ (%|2  % ] Re(w)a,oQe/aa,)J^1(k1a)lE|2[  eE^ 
dr|N^  \c  I yQ2  4irnj)inczYo  7 kmc^y^a 

generalizing  the  scaled  variables  to 

k2(3e/3u))J2+1(k1a)|E|2 

U 4irn,  my  3a 

b o 


kmc2y  a 
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mY3  (a)2y 

1 I! | 

2 fi.  N 


results  in 


nr  [q  + ui  = £ 

37  [K  + 2Re(w)u]  = £q 


which  are  the  analog  of  Eqs.(32)  and  (33)  and  in  the  limit  y - 1 reduce 

o 

identically  to  them. 

Next  we  find  the  connection  between  the  momentum  spread  A and  K,  where 

2 I 0VP)2 

*2  5 i — — <s 


or  in  scaled  variables 
r P2a2 

*2  - \ s?jr 

resulting  in 


P ay  * , 

* - - V- 

o 


hence 


S = i (2K-q2)^ 


To  check  the  relativistic  formulation  we  have  compared  its  predictions 
with  the  computer  study12  of  Lampe  and  Sprangle  in  which  the  full  relativistic 
dynamics  of  the  beam  particles  is  followed  for  the  case  Eq  = 0,  in  analogy 
with  the  nonrelativistic  work  of  O'Neil,  et  al.8  One  finds  again  that  the 
present  formulation  yields  results  which  average  over  the  relativistic  trapped 
particle  oscillations,  as  in  Fig.  1.  In  addition,  one  uncovers  a nontrivial 
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dependence  of  the  saturation  amplitude  on  the  value  of  the  parameter  s. 

The  existence  of  such  a dependence  was  originally  pointed  out  in  a computer 
simulation13  by  Thode  and  Sudan,  and  subsequently  investigated12  by  Lampe 
and  Sprangle  for  the  single  wave  problem,  and  by  Schamel,  et  al.17  in  the 
context  of  parametric  excitation  of  ion  density  fluctuations  by  the  single 
wave. 

The  solid  curve  in  Fig.  14  exhibits  the  predicted  dependence  of  the 

scaled  saturation  level  (s/4)u  on  s.  It  is  seen  that  for  small  s the  curve 

M 

increases  linearly  with  s,  as  expected,  and  for  s > 1.0  it  follows  a much 
slower  dependence,  as  found  earlier  by  Lampe  and  Sprangle  whose  results  are 
indicated  by  the  circles  in  Fig.  14.  It  is  observed  in  Fig.  14  that  the 
spatially  averaged  formulation  yields  an  answer  which  agrees  quantitatively 
with  the  results  of  the  exact  dynamical  calculation,  but  averages  over  the 
small  oscillations  present  in  the  more  detailed  description. 

After  establishing  the  predictive  power  of  the  formulation  we  have  pro- 
ceeded to  investigate  the  existence  of  the  clamping  behavior  in  the  relati- 
vistic regime.  As  is  illustrated  in  Fig.  15  for  £ = 0.1,  u(0)  = 1.0  * 10"6, 
s = 0.5,  one  finds  that  a runaway  relativistic  beam  can  experience  clamping, 
but  in  general  the  time  dependence  of  the  phenomena  exhibits  more  abrupt 
changes  than  in  the  nonrelativistic  case,  e.g.,  Fig.  9.  The  effects  of  multi- 
mode  interactions  and  sideband  generation  have  also  been  included  in  the 
formulation,  and  their  behavior  has  been  investigated.  However,  the  results 
do  not  differ  significantly  from  those  discussed  in  Sec.  IV  for  the  nonrelatis- 
tic  problem,  consequently  they  need  not  be  presented  here. 


and  p(x)  = p where  pQ  must  be  obtained  self-consistently  from  the  disper- 
sion relation. 

Equation  (65)  describes  the  existence  of  a clamped  state  in  which  the 
momentum  push  of  the  external  electric  field  is  transferred  to  the  cavity 
mode  and  subsequently  this  momentum  is  dissipated  through  the  appropriate 
damping  mechanisms  acting  on  the  mode  (e.g.,  wall  damping,  Landau  damping, 
radiation  leakage).  After  finding  the  existence  of  such  a dynamical  equi- 
librium one  must  inquire  into  its  stability.  The  stability  properties  are 
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determined  from  the  evolution  of  small  perturbations  p',  u' » such  that 
P * P0  ♦ P'»  u = uq  + u',  p'<<  P0  and  u'<<  uq.  The  corresponding  linearized 
system  is 

37  P'  - - 2r0u'  * 2 V’  (66) 

o 


2(r0-D)u'  . 2 (-g)  u0p' 
r o 


(67) 


seeking  a solution  p'  = p' (O)exp(yi) , u'  = u'(0)exp(yt)  results  in 


- fSKw  rJ 1 


- 4D(f)o“„/  <“) 


From  the  dependence  of  r on  p illustrated  in  Fig.  3 one  observes  that  there 
are  two  possible  values  of  pQ  such  that  Tq  < D,  as  required  by  the  steady 
state  solution.  One  root  corresponds  to  pQ  < 0 having  (3r/3p)Q  > 0,  while 
the  other  has  p > 0 and  (3T/3p)  < 0.  The  p < 0 root  is  stable  and  can 
support  damped  oscillations  about  the  stable  point  if 


2 (!£  (£♦  o)>  t(S  (4^)  + ^ 


3p 


Zp'  v 2D 
v o 


£ +oJ 


(69) 


The  pQ  > 0 root  is  stable  provided  that 


D > (.£*0)  [|(~)  \/2a]h  (70) 

v o 

For  large  values  of  the  external  electric  field  £ the  excitation  of 
the  cavity  mode  is  not  able  to  clamp  the  beam,  as  shown  earlier  in  Fig.  8. 

In  this  case  the  beam  momentum  p increases  monotonically  in  time,  hence 
the  condition  r > D is  eventually  attained  and  the  amplitude  of  the  cavity 
mode  grows.  The  time  evolution  of  this  process  can  be  described  analytically 
by  approximating  r(p)  s X | p | where  A is  a constant  that  can  be  determined 
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from  the  asymptotic  behavior  of  the  exact  T shown  in  Fig.  3.  In  this  regime 
P(t)  = - P(0)  + to  lowest  order,  thus  the  time  evolution  of  u is 

approximately  determined  by 


4-  u + 2 < D r 

| [|P(0)  |-£t]  '2 

in  the  region  p(t)  < 0. 

The  solution  of  Eq.  (71)  takes  the  form 


u = o 


(71) 


u(t)  = u(0) (exp  [-/QdT ' g (t ’ ) ] } { 1 + 2D/odT'exp[/odT"g(T")] } (72) 


with 


g(T) 


o.J 


| [Ip(0)|-£t] 

After  performing  the  required  integrations  one  obtains 

u(t)  = u(0)exp{-2DT-(4X/<f)[(|p(T)|)!s-(|p(0)|)li]}  (g) 

{1  + F(T)exp{(2X2/«Tn)  ♦ [2D|p(0)|/<f  ] - (4A/£  ) ( |p  (0)  | )**} } 


(73) 


(74) 


with 


F(x)  = A ( 2tt/ £ {erf[y(t)]  - erf[y(0)]} 


+ exp [-  y2(t)]  - exp [-  y2 (0) ] 


(75) 


and 


y(T)  = (2D/<f’)!s(|p(T)|)!s  - A (2/D<f*)^ 


(76) 


where  erf  refers  to  the  error  function. 


The  behavior  predicted  by  Eq.  (74)  is  illustrated  in  Fig.  16  for 
p(0)  = - 20.0,  £ = 1.0,  A = 0.87,  D = 0.25.  It  is  seen  in  Fig.  16  that 


I 


as  the  beam  momentum  approaches  the  threshold  point  (i.e.,  D - r = 0)  the 
amplitude  of  the  cavity  mode  increases  linearly  in  time;  beyond  the  thres- 
hold the  amplitude  grows  exponentially,  as  expected. 


-35- 


VII.  CONCLUSIONS 

This  study  illustrates  some  of  the  basic  processes  associated  with  the 
interaction  of  a runaway  electron  beam  with  collective  modes,  thus  broadening 
the  present  understanding  of  runaway  behavior  as  well  as  of  the  beam-plasma 
interaction.  The  principal  effect  found  is  that  a significant  amount  of  the 
push  provided  by  the  external  DC  electric  field  can  be  transformed  through 
the  beam  nonlinearities  to  the  form  of  wave  momentum,  hence  giving  rise  to 
the  clamping  of  the  beam.  This  phenomenon  has  two  consequences : (1)  it 

limits  the  maximum  energy  of  runaways;  (2)  it  provides  a mechanism  for 
converting  DC  to  AC  energy. 

In  solving  the  model  problem  a spatially  averaged  description  has  been 
introduced  which  permits  the  description  of  complicated  events  which  would 
normally  require  the  use  of  computer  simulation.  When  the  spatially  averaged 
technique  is  applied  to  the  pure  beam-plasma  interaction  it  is  found  that  the 
results  are  in  excellent  agreement  with  previous  computer  studies,  thus  it 
provides  additional  insight  into  the  basic  physics  behind  those  studies. 
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FIGURE  CAPTIONS 

Time  evolution  of  the  scaled  wave  amplitude  4T  and  scaled 
average  momentum  p for  a single  mode  and  Eq  = 0.  Solid 
curves  are  the  results  of  the  present  calculation,  dashed 
curve  is  the  result  of  O'Neil,  et  al.  (Ref.  8). 

Time  evolution  of  the  scaled  wave  amplitude  Su  and  scaled 
average  momentum  p for  a single  mode  and  Eq  = 0,  but  setting 
S = 0 explicitly.  To  be  compared  with  Fig.  1. 

Dependence  of  the  scaled  growth  rate  r and  scaled  frequency 
shift  0 on  p for  a cold  beam  S = 0,  and  for  S = 0.5. 

Self-consistent  time  evolution  of  the  scaled  growth  rate  r, 
scaled  frequency  shift  0,  average  kinetic  energy  K,  and 
thermal  spread  S for  the  behavior  shown  in  Fig.  1. 

Time  evolution  of  sideband  decay  for  Eq  = 0,  k^/k^  = 0.7, 
and  for  62  = 2.0,  1.5,  and  0.5. 

Time  evolution  of  the  scaled  growth  rate  r and  frequency 
shift  D during  the  sideband  decay  for  Eq  = 0,  k2/kj  = 

62  = 0.5.  To  be  compared  with  Fig.  5. 

Cascading  of  the  spectrum  from  a high  phase  velocity  mode  k^ 
to  a lower  phase  velocity  mode  k^  and  corresponding  evolution 
of  the  average  momentum  p and  thermal  spread  S.  Uj(0)  = u2(0) 

= 1.0  x lO-2,  p(0)  = 2.0,  Eq  = 0,  k2/k1  = 0.7,  62  = 3.0. 

Time  evolution  of  the  wave  energy  u^  and  average  momentum  p 
for  a runaway  case  £ = 1.0;  single  mode  with  u^(0)  = 1.0  x 10~3. 


FIG.  8 
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FIG.  9 


FIG.  10 


FIG.  11 


FIG.  12 


FIG.  13 


FIG.  14 


FIG.  15 


FIG.  16 


Example  of  clamping  behavior.  Time  evolution  of  wave  energy 
Uj  and  average  momentum  p for  £ = 1.0;  single  mode  with 
u^(0)  = 5.0  x 10‘3.  To  be  compared  with  Fig.  8. 

Time  evolution  of  the  effective  force  dp/dr  during  clamping 
and  self-consistent  behavior  of  the  kinetic  energy  K and 
thermal  spread  S for  the  case  shown  in  Fig.  9. 

Example  of  clamping  behavior  following  an  overshoot  of  the 
cavity  mode,  u^  is  the  wave  energy  and  p the  average  momentum. 
£=  1.0,  p(0)  = - 2.5,  S(0)  = 0,  u1  (0)  = 5.0  x 10'3. 

Example  of  clamping  behavior  in  the  presence  of  two  modes. 
k2/kj  = 0.7,  62  = 3.0,  £ = 1.5,  u^O)  = u2(0)  = 1.0  x 10'2, 

S(0)  = 0. 

Excitation  of  two  cavity  modes  by  a runaway  warm  beam.  £ = 1.5, 
u^O)  = u2(0)  = 5.0  x 10-2,  k2/k1  = 0.7,  «2  = 3.0,  S(0)  = 1.0. 

Dependence  of  scaled  saturation  energy  si^/4  on  the  parameter 
s for  the  relativistic  single  mode  Eq  = 0 problem.  Solid 
curve  is  the  prediction  of  the  present  model,  the  circles  are 
the  computer  simulation  results  of  Lampe  and  Sprangle  (Ref.  12). 

Example  of  clamping  behavior  in  the  relativistic  regime  for  a 


single  mode  case 


. <5  = 0.1, 


s = 0.5. 


Time  evolution  of  the  wave  energy  u and  total  growth  rate  D-r 
for  a damped  cavity  mode,  as  predicted  by  perturbation  theory. 
D = 0.25, £ = 1.0,  p(0)  = - 20.0,  X = 0.87. 
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